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ABSTRACT: In the work the matrix methods of flood risk analysis and assessment were presented. 
Also the issue of flood risk in terms of its regulations was presented. The flood risk analysis and 
assessment were performed with the use of risk matrix for one of the Subcarpathian commune. 
For this purpose, flood risk maps were prepared and potential flood losses were calculated in accord- 
ance with existing legislation. 
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Introduction 


Flood due to possible negative effects should be included in crisis situa- 
tion understood in accordance with art. 3 of the Act of April 26, 2007 on Cri- 
sis Management (Journal of Laws 2007 No. 89 item 590) as "a situation 
affecting negatively the level of safety of people, property of considerable size 
or the environment, causing significant restrictions on the operation of rele- 
vant public administration bodies due to the inadequacy of the forces and 
possessed means” (Dz.U. 2002 nr 62 poz. 558;Dz.U. 2017 poz. 209). 

The Floods Directive contains important issues related to the adopted 
flood risk strategies. In order to create the basis of these strategies, it is nec- 
essary to develop flood risk maps and, on their basis, to develop flood risk 
management plans (Ahmad et al., 2013; Blažkova, Beven, 2002; Brocca, 2013; 
Flores-Montoya et al., 2016; Róthlisberger et al., 2017). In the latter, Member 
States should refrain from taking measures and implementation of measures 
that could increase the risk of flooding in the other Member States. And even, 
in accordance with the principle of solidarity, the decision should be made 
jointly and act along the entire course of the river (Graniczny, Mizerski, 2007). 

For each river basin district, Member States should prepare a prelimi- 
nary flood risk assessment, containing at least: river basin maps, description 
of floods that occurred in the past if the likelihood of similar floods occur- 
rence in future is big, description of floods from the past if similar events in 
the future can have negative consequences. On the basis of an initial flood 
risk assessment, Member States identify areas where there is a high risk of 
flooding or the occurrence of such a risk is likely. For each river basin district 
flood hazard maps and flood risk maps should be prepared. Flood hazard 
maps show geographic areas where the probability of flooding is low or 
extreme events are likely to occur. For each of the above scenarios, the follow- 
ing elements are depicted on the maps: flood range, water depth and, depend- 
ing on the needs, the water flow rate. Flood risk maps present possible nega- 
tive consequences of flooding and include: estimated number of inhabitants 
affected by floods, type of economic activity in the area affected by flooding 
and other useful information depending on the Member State. On the basis of 
flood risk maps Member States develop flood risk management plans and set 
objectives for flood risk management “with particular emphasis on reducing 
the potential negative consequences of flooding for human health, the envi- 
ronment, cultural heritage and economic activity” (Bartnik, Jokiel, 2012; 
Dyrektywa 2007/60/WE). 

From the known methods of risk analysis and assessment to perform 
a flood risk analysis, the most appropriate is risk matrix that is based on 
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a basic risk definition (Cuellar, McKinney, 2017; Haimes, 2009; Lubowiecka, 
Wieczysty, 2000; Rak, 2004, 2008, 2009, 2010, 2014; Rincón etal., 2018; Sieg 
etal., 2017; Tchórzewska-Cieślak, 2008; Zhou et al., 2018). In order to create 
flood risk maps, the Geographic Information Systems (GIS) should be applied 
to the numerical terrain model (NMT) of water table elevations during floods 
obtained as a result of mathematical hydrological modelling (Tokarczyk et al., 
2012). To set an example to show the issue, the exemplary matrices of flood 
risk for the commune X located in the Subcarpathian province, were pre- 
sented. As a determinant of losses during floods, the number of inhabitants 
residing in the area at risk of flooding according to the flood hazard map was 
adopted. In this way, flood risk matrices were prepared for X commune, tak- 
ing into account the number of people at risk of flooding in relation to the 
flooding probability. 

The Water Law defines the flood risk as “the combination of the probabil- 
ity of a flood event and of the potential adverse consequences for human 
health, the environment, cultural heritage and economic activity associated 
with a flood event” (Dz.U. 2017 poz. 1566). 

Flood protection is the task of government and self-government adminis- 
tration bodies assisted by water users. The protection against floods is per- 
formed taking into account flood hazard maps (Armenakis et al., 2017; Daw- 
son et al., 2011; Haberlandt et al., 2008; Ouma, Tateishi, 2014; Simonovic, 
2009; Zischg, 2018), flood risk maps and flood risk management plans, 
through e.g. multi-index evaluation for flood disaster (Dou et al., 2018), 
building ANN-Based Regional Multi-Step-Ahead Flood models (Chang et al., 
2018), using a response curve approach (Murdock et al., 2018), sketch maps 
(Klonner et al., 2018), cloud-model-based method for risk assessment (Yang 
et al., 2018), application of a bayesian approach to dynamic assessment of 
flood (Wu et al., 2018), a parametric distance function approach (Zheng et al., 
2018). 

For their preparation, it is necessary to conduct preliminary flood risk 
assessment, taking into account the method of calculating the value of poten- 
tial flood losses in individual classes of land use, for the purposes of develop- 
ing flood risk maps (Dz.U. 2013 poz. 104). The paper proposes the use of 
matrix methods for the analysis and assessment of flood risk in accordance 
with the applicable legislation. 
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Research methods 


Matrix methods for risk analysis 


The risk matrix is a combination of a point scale probability of the unde- 
sirable event with a scale of consequences. The two-parameter risk estima- 
tion matrix is one of the simplest. Point values defining the risk categories 
were obtained using the dependence (Rak, 2004): 


r=P.C (1) 
where: 
P - the probability of the undesirable event, 
C - consequences of this event. 


Considering the complexity of the systems, a three-parameter risk matrix 
was proposed. Additional parameters included in the matrix are: exposure to 
hazard - E or vulnerability - V. According to the above parameters, the 
numerical risk assessment is their product (Rak, 2004; Rak, Tchórzews- 
ka-Cieślak, 2005): 


r=P-C-E (2) 
or 
r=P-C-V (3) 


Along with the development of modern techniques, almost all the areas 
of life are equipped with various types of security and monitoring systems, in 
order to increase their safety and reliability. Therefore, it seems advisable to 
include in the matrix a fourth parameter defining the amount of protection. 

A four-parameter matrix for risk estimation based on the formula is pro- 
posed (Rak, 2004; Rak, Tchórzewska-Cieślak, 2005): 

P-C-V 


GE (4) 


where: 
O - protection against threats. 


In addition to protection in the risk matrix, the risk exposure (E) known 
from the three-parameter matrix can be taken into account, in that way we 
obtain a five-parameter risk estimation matrix along with the formula (Rak, 
2004; Rak, Tchórzewska-Cieślak, 2005): 

P-C-M-E 


Pag ©) 


where: 
O - protection point weight. 
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Each time, for individual parameters, the size level is assigned according 
to the adopted point scale, e.g. low - N = 1, medium - S = 2, high - W=3. 

In this way, a punctual scale of risk measures is obtained in numerical 
form, which is the basis for risk assessment. 

The acceptable level of risk is determined by introducing criterion values 
for each level of tolerable, controlled and unacceptable risk. 


Implementation of risk matrices in the assessment of flood risk 


Adaptation of a two-parameter matrix 


On the basis of the number of endangered inhabitants, consequences cat- 
egories (C) presented in table 1 were adopted. In table 2 the categories of 
probability (P) of flood occurrence are presented (Zygmunt, 2017). 


Table 1. Categories of consequences — C 


Number of endangered inhabitants Point weight Description 
0 1 No danger 
1-10 2 Negligible 
11-100 3 Marginal 
101-200 4 Significant 
>200 5 Serious 


Source: authors own work. 


Table 2. Probability categories — P 


Probability of flood occurrence Point weight Description 
0,5% 1 Improbable 
1% 2 Unlikely 

5% 3 Occasional 
10% 4 Possible 
50% 5 Frequent 


Source: authors own work. 


After combining the probability and consequences categories, a risk 
matrix was obtained according to the dependence 1, which was presented in 
the form of table 3. 
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Table 3. Risk matrix 


PROBABILITY - P 


0,5% % 5% 10% 50% 
POINT WEIGHT 
Improbable Unlikely Occasional Possible Frequent 
3 4 5 
No threat 1 1 2 3 4 5 
© 
«, Negligible 2 2 4 6 8 10 
Lu 
z Marginal 3 3 6 9 12 15 
& Significant 4 4 8 12 16 20 
z 
S Serious 5 5 10 15 20 25 


Source: authors own work. 


Depending on the obtained results in the risk matrix theflood risk was 
divided into: 
e weight from 1 to 4 - tolerable risk, 
e weight from 6 to 10 - controlled risk, 
e weight from 12 to 25 - unacceptable risk. 


Adaptation of a three-parameter matrix 


Considering the vulnerability (V) related to flood risk zones, the risk 
assessment may be determined by multiplying the probability parameter 
(P), the consequences (C) and the vulnerability to the risk (V), according to 
the equation 3. 
The following scale and weight for individual parameters were proposed: 
e point weights for the probability parameter - P: 
- improbable event, the probability of flood occurrence 0,5%; with 
a weight 1, 

- unlikely event, the probability of flood occurrence 1%; with a weight 2, 

- occasional event, the probability of flood occurrence 5%; with 
a weight 3, 

- possible event, the probability of flood occurrence 10%; with a weight 
of 4, 

- frequent event, the probability of flood occurrence 50%; with a weight 
of 5. 

e point weights for the consequences parameter - C: 
- no threat, number of endangered inhabitants 0; with a weight of 1, 
- negligible, number of endangered inhabitants 1-10; with a weight of 2, 
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marginal, number of endangered inhabitants 11-100; with a weight 
of 3, 
significant, number of endangered inhabitants 101-200; with a weight 
of 4, 
serious, number of endangered inhabitants > 200; with a weight of5. 


e point weights for the vulnerability parameter — V (Nachlik et al., 2000; 
Ozga-Zielińska et al., 2003; Radczuk et al., 2007): 


flood threat zone (ZS) - this area is for the most part of the year con- 
stantly flooded with water up to the average of the highest flows in mul- 
ti-year (SWQ); for a water depth range of < 0,5 m, with a weight of 1, 
high flood hazard zone (ZW) - this is the area above the ZS zone and 
with the upper limit at the level of water from the maximum flow 
with the probability of exceeding 1% of the so-called hundred-year 
water (Qmax 1%); for a water depth range of 0,5 < h < 2 m, with 
a weight of 2, 

significant flood zone (ZZ) - this is the area above the ZW zone with 
the upper limit at the level of water caused by the flow with half the 
maximum reliable flood (0,5xMWW); for a water depth range 2 < h < 
4 m, with a weight of 3, 

small flood hazard zone (ZM) - this is the area above the ZZ zone up 
to the maximum reliable flood (MWW); for a water depth range 
> 4 m, with a weight of 4. 


In this way, we get the quantitative risk matrix presented in table 4. 


Table 4. Risk matrix 


PROBABILITY - P 


0,5% 1% 5% 10% 50% 
POINT WEIGH 
Improbable Unlikely Occasional Possible Frequent 
| + 52> 535 4-1. 022 UB A AA Bo Ah 290042 As ADP BA 
No threat 1 12 3 4 2 4 6 8 3 6 9 12 4 8 12 16 5 10 15 20 
Negligible 2 2 4 6 8 4 8 12 16 6 12 18 24 8 16 24 32 10 20 30 40 
z Marginal 3 3 6 9 12 6 12 18 24 9 18 27 36 12 24 36 48 15 30 45 60 
g Significant 4 4 8 12 16 8 16 24 32 12 24 36 48 16 32 48 64 20 40 60 80 
& Serious 5 5 10 15 20 10 20 30 40 15 30 45 60 20 40 60 80 25 50 75 100 


Source: authors own work. 


In this way, point scale measures of risk in the numerical range from 1 to 
100, for the following risk scale from 1 to 4 (tolerable risk), from 5 to 10 
(controlled risk), from 12 to 100 (unacceptable risk) were obtained. 
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Adaptation of a four-parameter matrix 


In order to assess the system of protection against threats, a four-param- 
eter risk assessment matrix, according to formula 3, was implemented. 

The description of the protection category is as follows (Ozga-Zielińska et 
al., 2003): 

e low level of protection, basic zone — ZS and ZW hazard zones - it is a zone 
permanently protected against flooding, because the flow corresponding 
to the upper limit of the ZW threat zone is designed for technical flood 
protection solutions, e.g. flood banks, point weight 1, 

e medium level of protection, extraordinary zone — ZZ threat zone — protec- 
tion in this zone is introduced in the situation of a threat in the area where 
the population feels safe as it is protected by solutions applied at the bor- 
der of the primary zone, point weight 2, 

e high degree of protection, extreme zone - danger zone ZM - protection in 
this zone is usually activated after occurrence of events causing extreme 
high flow, point weight 3. 

Other parameters are presented in point 3.2.The presented scale of point 
weights is a proposal for the initial risk assessment, in this way the following 
risk categories can be assumed: 

e tolerable, point scale from 0,33 to 3, 

e controlled, point scale from 3 to 20, 

e unacceptable, point scale from 20 to 100. 


Determination of potential flood losses 


The estimated number of residents likely to be affected by the flood is 
assumed to be the number of people registered in the area. Residential build- 
ings and objects of special social significance are distinguished between 
those “for which the water depth is less than or equal to 2 m, and objects for 
which the water depth is greater than 2 m” (Dz.U. 2013 poz. 104). 

The values of potential flood losses for individual classes of land use are 
determined by the formula (Dz.U. 2013 poz. 104): 


4 
SP; => Sp; -Ai fori=1...9 (6) 
j=1 


where: 
Sp; - means the total values of potential unit losses for a given class of land use, 
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Sp; - means the value of potential unit losses for a given class of land use and the 
water depth range, 
A; - means the area occupied by a given land-use class. 


However the values of potential unit losses in classes are expressed as 
(Dz.U. 2013 poz. 104): 


Sp, =w, -f(h,) (7) 
where: 
Sp; - means the value of potential unit losses for a given class of land use and the 
water depth range, 
W, - means the value of the property in a given class of land use, 


f(h,)- means the value of the loss function relating to the depth of water with the 
loss of value of the property in a given class of land use. 


The values of the loss function are selected for each class of land use accord- 
ing to the water depth at the limits of the intervals of these depths from the tables 
contained in the regulation (Dz.U. 2013 poz. 104). In turn, the value of assets in 
a given class of land use is selected from tables depending on the voivodship. 

In table 5 the values of the loss function for residential areas which, 
according to the regulation, should be taken into account in order to link the 
water depth during floods with the loss of property value in a given class of 
land use, were presented. 


Table 5. The value of the loss function for residential areas. 


Water depth h [m] Value of the loss function f(h) [%] 
<0,5 20 
05<h<2 35 
2<h<4 60 
>4 95 


Source: authors own work based on (Dz.U. 2013 poz. 104). 


Application example 


The analysis was conducted on the basis of materials received from the office 
of the distinguished commune located in the Subcarpathian province and occu- 
pying an area of 200 km. For analysis, 12 villages of the considered commune 
were taken. The basic data were flood area coverage maps with a given probabil- 
ity of flooding. Using the ArcGIS Explorer version, a flood risk map was drawn up. 
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Table 6. Flood risk in the distinguished commune 


Flood threat 
Location number The probability The number of ; Number of 
of residents Percentage ratio 
of flood vulnerable households endangered people 
A 3563 due to the lack of a large watercourse flowing through the distinguished location, there are no people at risk of 
being exposed to a flood 
B 1399 despite the flow of the stream through the settlement, there are no people at risk of flooding 
C 896 no watercourse threatening the inhabitants, so there are no people at risk of flooding 
D 1415 due to the location of the buildings directly on the banks of the river and near the estuary, the number of 
people at risk of flooding was determined depending on the flood probability 
50% 0,3% 4 
0% 3 % 13 
5% 0 3% 45 
% 42 3% 188 
0,5% 57 8% 255 
E 4191 the stream flowing through the village is a threat to the inhabitants only in the case of extremely high flood 
0,5% 3 0% 14 
F 788 due to the lack of a large watercourse flowing through the village, there are no people at risk of a flood 
318 no watercourse and endangered inhabitants 
1355 the location of the village on the banks of the river poses a great threat to the residents, therefore the number 
of people at risk was determined 
50% 2 1% 0 
0% 20 1% 98 
5% 32 12% 57 
% 61 22% 298 
0,5% 91 33% 445 
50% 2 1% 0 
I 972 due to the lack of a watercourse in the village, there are no people exposed to the flood 
J 293 due to the proximity of the river, the number of endangered inhabitants is presented in the next lines 
0% 5 8% 22 
5% 13 20% 58 
% 30 45% 33 
0,5% 32 48% 42 
K 1749 there is a threat due to the proximity of the river 
0% 5 „3% 23 
5% 6 2% 27 
% 60 6% 273 
0,5% 70 8% 318 
L 3483 as in the case of the village of E, despite the flow of the stream, the flood risk is not too high and the number 


of endangered inhabitants is presented in the following 


lines 


5% 1 0,1% 5 
1% 4 % 19 
0,5% 8 % 39 


Source: author's own work. 
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On the basis of the flood hazard maps, the flood hazard number obtained 
in each location and the number of households, the percentage of all house- 
holds in a given location in danger of flooding, was calculated (table 6) (Zyg- 
munt, 2017). 

For particular location of the commune scales of the probability (P) and 
the consequences (C) were assumed and the value of flood risk (r) was calcu- 
lated, according to the formula 1, using the guidelines of the two-parameter 
matrix. Summary of the results together with the division of flood risk is pre- 
sented in table 7. 


Table 7. Flood risk in the concerned commune 


Location P C r Description 
A 1 1 Tolerable risk 
B 1 1 Tolerable risk 
C 1 1 Tolerable risk 
D 2 4 8 Controlled risk 
E 2 2 olerable risk 
F 1 olerable risk 
G 1 olerable risk 
H 2 5 0 Controlled ris 
1 olerable risk 
J 2 4 8 Controlled ris 
K 2 5 0 Controlled ris 
L 2 3 6 Controlled ris 


Source: authors own work. 


On the basis of the obtained results for the flood risk in particular locali- 
ties of the analysed commune, the estimated flood losses in the housing areas 
of the commune were calculated. 

The formula presented in the Regulation of the Minister of Environment, 
Minister of Transport, Construction and Maritime Economy, Minister of 
Administration and Digitization and Minister of the Interior of 21 December 
2012 on the development of flood hazard maps and flood risk maps was used 
for this purpose (Dz.U. 2017 poz. 1566). 

According to formula 5, the values of potential flood losses for residential 
areas and water depth ranges were calculated using formula 6. 
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The value of assets for residential housing areas, which according to the 
regulation is adopted for the calculation of flood losses depending on the 
province, for the Subcarpathian province is 201,25 PLN/m? (Dz.U. 2017 poz. 
1566). The obtained values of potential unit losses for residential areas and 
water depth ranges are presented in table 8. 


Table 8. Values of potential flood losses 


Water depth h [m] The value of the loss function f(h) [%] Sp; 
<0,5 20 40,25 
05<h<2 35 10,44 
2<h<4 60 120,75 
>4 95 191,19 


Source: author's own work based on (Dz.U. 2017 poz. 1566). 


In order to estimate the area occupied by residential buildings in individ- 
ual location of the considered commune flooded with water during the flood, 
with the probability assumed in the flood risk matrix with the division into 
water depth ranges, according to table 7, flood hazard maps were used. Anal- 
ogously to the earlier proceedings, the map was enlarged to a size that makes 
it possible to distinguish individual households and the number of flooded 
householdsin the given depth ranges was calculated. Places in which tolera- 
ble flood risk occurred have not been taken into account in further calcula- 
tions. Therefore flood losses were calculated for five locations in the analysed 
commune, in which the controlled flood risk was found (table 9). 

For individual towns where flood losses occurred, the weights were cal- 
culated for the flood risk value (r), according to the formula 4, using a four-pa- 
rameter matrix. The summary of the results together with the division of 
flood risk is presented in table 10. 
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Table 9. Number of households in water depth ranges 


Location Ranges of water depth [m] 
<0,5 05<h<2 2<h<4 >4 
Number of households in water depth ranges 21 21 0 0 
D Percentage of the water depth ranges [%] 6,65 6,65 0,00 0, 
The value of flood losses in thousand [PLN] 1540 2695 0,00 0, 
Number of households in water depth ranges 35 22 4 0 
H Percentage of the water depth ranges [%] 12,64 7,94 1,44 0, 
The value of flood losses in thousand [PLN] 3472 3819 1190 0, 
Number of households in water depth ranges 21 9 0 0 
J Percentage of the water depth ranges [%] 31,82 13,64 0,00 0, 
The value of flood losses in thousand [PLN] 194 2 494 1871 0, 
Number of households in water depth ranges 25 32 3 0 
K Percentage of the water depth ranges [%] 6,49 8,31 0,78 0, 
The value of flood losses in thousand [PLN] 2328 5215 838 0, 
Number of households in water depth ranges 4 0 0 0 
L Percentage of the water depth ranges [%] 0,56 0,00 0,00 0, 
The value of flood losses in thousand [PLN] 274 0,00 0,00 0, 
Source: author's own work based on (Dz.U. 2017 poz. 1566) 
Table 10. Flood risk in the concerned commune 
Location P C V r Description 
D 2 4 1 8 Controlled risk 
2 4 2 16 Controlled risk 
H 2 5 1 10 Controlled risk 
2 4 2 8 Controlled risk 
2 4 3 12 Controlled risk 
J 2 4 1 8 Controlled risk 
2 3 2 12 Controlled risk 
K 2 5 1 10 Controlled risk 
2 5 2 20 Controlled risk 
2 3 3 18 Controlled risk 
L 2 3 1 6 Controlled risk 


Source: authors own work. 
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The estimated potential flood losses for residential areas in the con- 
cerned municipality at the occurrence of a flood with a probability of 1%, i.e. 
the floods of the century, are over PLN 25 million. Places in the order of the 
highest losses are: H, K, J, D and L. 


Conclusions 


Flood risk analysis and assessment aims to protect people, property and 
the environment from the effects of flood, it is developed in the form of flood 
hazard maps, flood risk and flood risk management plans. 

The following conclusions were proposed after performed flood risk 
analysis: 

e on the basis of the analysis of table 9, it was found that the flood risk in 
the considered commune depends on the distance of individual towns 
from the largest river flowing through the commune, 

e the unacceptable risk does not occur in any location in the concerned 
commune, 

e the controlled risk occurs in five locations of the commune: D, H, J, K and 
L. All these towns are in the vicinity of the river, 

e the tolerable risk occurs in the other seven locations of the commune: 
A, B, C, E, F, G and I. None of these locations are adjacent to the river, 

e the concerned commune is a relatively safe place in terms of flood risk 
when taking into account inhabited areas. For areas used for agriculture, 
flooding is likely to be more severe, because in these areas water depth 
will be larger, 

e the largest flood hazard in the concerned commune is found in the loca- 
tion near the only bridge in the commune. In case of “the flood of the 
century” nearly half of the households in the J location will be flooded 
and almost 1/3 of the households in the neighbouring location H will be 
flooded, 

e aslightly smaller threat is found in the village of K, also located near the 
bridge and in the village of D in the vicinity of the estuary, because less 
than 1/5 of the farms will be flooded with hundred-year-old wate, 

e flood losses occur adequately to the flood risk: the higher the risk, the 
greater the losses. 

For the location, where the controlled flood risk was found, exemplary 
flood losses for residential areas were calculated in accordance with the 
guidelines included in the regulation on the preparation of flood risk maps 
and flood risk maps. The flood losses expressed in monetary values in the 
event of “the flood of the century” will amount to approx. 25 740 thousand 
PLN. 
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